A Search for Electron Neutrino Appearance at the Am 2 ~ 1 eV 2 Scale 
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The MiniBooNE Collaboration reports first results of a search for v e appearance in a beam. 
With two largely independent analyses, we observe no significant excess of events above background 
for reconstructed neutrino energies above 475 MeV. The data are consistent with no oscillations 
within a two-neutrino appearance-only oscillation model. 



This Letter reports the initial results from a search for 
Vy, — ► v e oscillations by the MiniBooNE Collaboration. 
MiniBooNE was motivated by the result from the Liquid 
Scintillator Neutrino Detector (LSND) experiment 



which has presented evidence for 



oscillations at 



the Ar 



1 eV scale. Although the Karlsruhe Ruther- 



ford Medium Energy Neutrino Experiment (KARMEN) 
observed no evidence for neutrino oscillations Q , a joint 
analysis Q showed compatibility at 64% CL. Evidence 
for neutrino oscillations also comes from solar-neutrino 
0, @, @, 0, S] and reactor-antineutrino experiments @, 
which have observed v e disappearance at Am 2 ~ 8 x 10 -5 
eV 2 , and atmospheric- neutrino 13, U, 13, r S| an d long- 
baseline accelerator- neutrino experiments [14|, [l5[ , which 
have observed disappearance at Am 2 ~3x 10 -3 eV 2 . 

If all three phenomena are caused by neutrino oscilla- 
tions, these three Am 2 scales cannot be accommodated 
in an extension of the Standard Model that allows only 
three neutrino mass eigenstates. An explanation of all 
three mass scales with neutrino oscillations requires the 



addition of one or more sterile neutrinos |16I| or further 
extensions of the Standard Model (e.g., [171]). 

The analysis of the MiniBooNE neutrino data pre- 
sented here is performed within a two neutrino 
appearance-only — > v e oscillation model which uses 
U/j, events to constrain the predicted v e rate. Other than 
oscillations between these two species, we assume no ef- 
fects beyond the Standard Model. 

The experiment uses the Fcrmilab Booster neutrino 
beam, which is produced from 8 GeV protons incident 
on a 71-cm-long by 1-cm-diamctcr beryllium target. The 
proton beam typically has 4 x 10 12 protons per ~ 1.6 
/Lts beam spill at a rate of 4 Hz. The number of pro- 
tons on target per spill is measured by two toroids in the 
beamline. The target is located inside a focusing horn, 
which produces a toroidal magnetic field that is pulsed 
in time with the beam at a peak current of 174 kA. Posi- 
tively charged pions and kaons, focused by the horn, pass 
through a 60-cm-diamctcr collimator and can decay in a 
50- m- long tunnel, which is 91 cm in radius and filled with 



2 



air at atmospheric pressure. 

The center of the detector is 541 m from the front 
of the beryllium target and 1.9 m above the center of 
the neutrino beam. There is about 3 m of dirt over- 
burden above the detector, which is a spherical tank of 
inner radius 610 cm filled with 800 tons of pure mineral 
oil (CH2) with a density of 0.86 g/cm 3 and an index of 
refraction of 1.47. The light attenuation length in the 
mineral oil increases with wavelength from a few cm at 
280 nm to over 20 m at 400 nm. Charged particles pass- 
ing through the oil can emit both directional Chcrcnkov 
light and isotropic scintillation light. An optical bar- 
rier separates the detector into two regions, an inner vol- 
ume with a radius of 575 cm and an outer volume 35 cm 
thick. The optical barrier supports 1280 equally-spaced 
inward-facing 8-inch photomultiplicr tubes (PMTs), pro- 
viding 10% photocathode coverage. An additional 240 
tubes are mounted in the outer volume, which acts as 
a veto shield, detecting particles entering or leaving the 
detector. Two types of PMT are used: 1198 Hamamatsu 
model R1408 with 9 stages and 322 Hamamatsu model 
R5912 with 10 stages. Approximately 98% of the PMTs 
have worked well throughout the data taking period. 

The experiment triggers on every beam spill, with all 
PMT hits recorded for a 19.2 /js window beginning 4.4 /js 
before the spill. Other triggers include a random trig- 
ger for beam-unrelated measurements, a laser-calibration 
trigger, cosmic-muon triggers, and a trigger to record 
neutrino- induced events from the nearby Neutrinos at the 
Main Injector (NuMI) beamline 11811 . The detector elec- 
tronics, refurbished from LSND [l9j|, digitize the times 
and integrated charges of PMT hits. PMT hit thresholds 
are ~0.1 photoclcctrons (PE); the singlc-PE time reso- 
lutions achieved by this system are ~ 1.7 ns and ~ 1.2 
ns for the two types of PMTs. One PE corresponds to 
~0.2 MeV of electron energy. Laser calibration, consist- 
ing of optical fibers that run from the laser to dispersion 
flasks inside the tank, is run continuously at 3.33 Hz to 
determine PMT gains and time offsets. Averaged over 
the entire run, the beam-on livetime of the experiment is 
greater than 98%. 

The energy spectrum peaks at 700 MeV and extends 
to approximately 3000 MeV. Integrated over the neutrino 
flux, interactions in MiniBooNE arc mostly charged- 
current quasi-elastic (CCQE) scattering (39%), neutral- 
current (NC) elastic scattering (16%), charged-current 
(CC) single pion production (29%), and NC single pion 
production (12%). Multi-pion and deep-inelastic scatter- 
ing contributions are < 5%. NC elastic scattering, with 
only a recoil nucleon and a neutrino in the final state, 
typically produces relatively little light in the detector 
and contributes only 3 events to the final background 
estimate. 

Tabic U shows the estimated number of events with 
reconstructed neutrino energy, E$ E , between 475 MeV 
and 1250 MeV after the complete event selection from 



TABLE I: The estimated number of events with systematic 
error in the 475 < E$ E < 1250 MeV energy range from all of 
the significant backgrounds, together with the estimated num- 
ber of signal events for 0.26% — > v £ transmutation, after 
the complete event selection. 



Process 


Number of Events 


Vy. CCQE 


10 ± 2 


f M e — » v^e 


7 ± 2 


Miscellaneous Events 


13 ±5 


NC 7T° 


62 ± 10 


NC A — > JV7 


20 ±4 


NC Coherent & Radiative 7 


< 1 


Dirt Events 


17 ±3 


v £ from [i Decay 


132 ± 10 


Ve, from K Decay 


71 ± 26 


v e from K% Decay 


23 ±7 


i/ e from 7r Decay 


3± 1 


Total Background 


358 ± 35 


0.26% Vy -> v e 


163 ±21 



all of the significant backgrounds, where E® E is deter- 
mined from the reconstructed lepton energy and angle 
with respect to the known neutrino direction. The back- 
ground estimate includes antineutrino events, which rep- 
resent < 2% of the total. Also shown is the estimated 
number of v e CCQE signal events for the LSND central 
expectation of 0.26% — > v e transmutation. Studies 
of random triggers have established that no significant 
backgrounds survive the analysis cuts other than those 
due to beam related neutrinos, which can be divided into 
cither ^-induced or j/ e -induced backgrounds. The small 
fraction of v e from /1, K, and tt decay in the beamline 
gives a background that is indistinguishable from oscil- 
lations except for the energy spectrum. CC events 
are distinguished from v e events by the distinct patterns 
of Cherenkov and scintillation light for muons and elec- 
trons, as well as by the observation of a delayed electron 
from the muon decay, which is observed > 80% of the 
time from CCQE events. NC ir° events with only a 
single electromagnetic shower reconstructed are the main 
^-induced background, followed by radiative A decays 
giving a single photon, and then neutrino interactions in 
the dirt surrounding the detector, which can mimic a sig- 
nal if a single photon, mostly from tt° decay, penetrates 
the veto and converts in the fiducial volume. 

We use PMT charge and time information in the 19.2 
jj,s window to reconstruct neutrino interactions and iden- 
tify the product particles. This time window is defined 
as an "event" and is divided into "subevents", collec- 
tions of PMT hits clustered in time within ~ 100 ns. A 
Vy CCQE event with a muon stopping within the tank 
may have two subevents: the first subevent from par- 
ticles produced at the neutrino interaction, the second 
from the muon decay to an electron. A v e CCQE event 
has a single subevent. 
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To ensure stable, well-targeted beam at full horn cur- 
rent, it is required that the two monitoring toroids agree 
to within 5%, the estimated transverse containment of 
the beam in the target be greater than 95%, and the 
measured horn current be within 3% of its nominal value. 
The event time at the detector must be consistent with 
the beam delivery time (both determined by GPS), and 
the event must pass a number of data integrity checks. 
The beam quality requirements reject 0.7% of the events, 
while the detector time and quality requirements remove 
a further 1.8%, with the remaining data corresponding 
to (5.58 ± 0.12) x 10 20 protons on target. 

Next, events with exactly one subevent (as expected 
for v e CCQE events) are selected. By requiring that the 
subevent have fewer than 6 hits in the veto and more 
than 200 hits in the main tank (above the muon-decay 
electron endpoint), entering cosmic-ray muons and their 
associated decay electrons are eliminated. The average 
time of hits in the subevent is required to be within the 
beam time window of 4-7 These cuts yield a cosmic 
ray rejection of greater than 1000:1. 

After these initial cuts, the surviving events are re- 
constructed under four hypotheses: a single electron-like 
Cherenkov ring, a single muon-like ring, two photon-like 
rings with unconstrained kinematics, and two photon- 



like rings with M, 



77 



n^o (see Fig. []}. Photon- like 
rings are assumed to be identical to electrons, but al- 
lowed to be independently displaced from the neutrino 
interaction vertex. The reconstruction uses a detailed 
model of extended-track light production and propaga- 
tion in the tank to predict the charge and time of hits on 
each PMT. Event parameters are varied to maximize the 
likelihood of the observed hits, yielding the vertex posi- 
tion and time of the event and the direction, energy, and, 
for photons, the conversion distance of the ring(s). For i/ e 
events, the event vertex, direction, and energy are recon- 
structed on average with resolutions of 22 cm, 2.8°, and 
11%, respectively, while NC tt° events are reconstructed 
with a 7T° mass resolution of 20 MeV/c 2 . 

The final analysis cuts were designed to isolate a sam- 
ple of ^-induced events that were primarily CCQE. The 
only data that were used in developing the analysis were 
samples that Monte Carlo (MC) simulation had indicated 
could not contain a significant number of — > v e oscil- 
lation events. We require that the electron-hypothesis 
event vertex and muon-hypothesis track endpoint occur 
at radii < 500 cm and < 488 cm, respectively, to en- 
sure good event reconstruction and efficiency for pos- 
sible muon decay electrons. We require visible energy 
E vis > 140 MeV. We then apply particle identification 
(PID) cuts to reject muon and 7r° events. These are E V i S - 
dependent cuts on \og(L e / L^), log(L e /L n o), and Af 77 , 
where L ei L^, and L^o are the likelihoods for each event 
maximized under the electron 1-ring, muon 1-ring, and 
fixed-mass 2-ring fits, and M 77 is from the unconstrained 
two-ring fit. These also enhance the fraction of CCQE 
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FIG. 1: Events in MiniBooNE are reconstructed as either a 
muon event, an electron event, or a tv° event. 



TABLE II: The observed number of u e CCQE candidate 
events and the efficiency for — > v £ CCQE oscillation events 
after each cut is applied sequentially. 



Selection 


#Events 


u e CCQE Efficiency 


Cosmic Ray Cuts 


109,590 


100% 


Fiducial Volume Cuts 


68,143 


55.2 ± 1.9% 


PID Cuts 


2037 


30.6 ± 1.4% 


475 < E® B < 1250 MeV 


380 


20.3 ± 0.9% 



events among the surviving electron candidates. Table HT1 
shows the observed number of v e CCQE candidate events 



and the efficiency for 



CCQE oscillation events 



after each cut is applied sequentially. A total of 380 data 
events remains after the complete selection. 

Detailed Monte Carlo simulations of the beam and de- 
tector were used to make initial estimates of the flux and 
detector efficiencies. The Booster neutrino beam flux at 
the detector is modeled using a GEANT4-bascd simu- 
lation [2(| of the beamline geometry. Pion and kaon 
production in the target is parametrized [2l| based on 
a global fit to proton-beryllium particle production data 
[22|. The kaon flux has been cross-checked with high- 
energy events above 1250 MeV and with an off-axis muon 
spectrometer that viewed the secondary beamline from 
an angle of 7°. This detector determined the flux of 
muons with high transverse momentum, which originate 
mostly from kaon decays, to be consistent with the MC 
predictions. 

The v3 NUANCE [23| event generator simulates neu- 
trino interactions in mineral oil. Modifications are made 
to NUANCE which include adjustment of the axial form 
factor of the nucleon for quasi-elastic scattering, the Pauli 
blocking model, and coherent pion production cross sec- 
tions based on fits to MiniBooNE data 24|. In addi- 
tion, the final state interaction model has been tuned to 
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reproduce external pion-carbon rescattcring data [25[ , an 
explicit model of nuclear de-excitation photon emission 
for carbon has been added [26| , and the angular correla- 
tions for A decay are modified to be in accord with the 
model of Rein and Sehgal j2?| ■ 

Particles from NUANCE-generated final states are 
propagated through a GEANT3-based simulation [28[ 
of the detector, with the subsequent decays, strong, 
and electroweak interactions in the detector medium in- 
cluded. Most particular to MiniBooNE, the emission of 
optical and near-UV photons via Cherenkov radiation 
and scintillation is simulated, with each photon individ- 
ually tracked, undergoing scattering, fluorescence, and 
reflection, until it is absorbed [2^]. Small-sample mea- 
surements of transmission, fluorescence, and scattering 
are used in the model. Muon decay electrons are used 
to calibrate both the light propagation in the detector 
and the energy scale. The amount of scintillation light 
is constrained from NC elastic scattering events. The 
charge and time response of the electronics is simulated, 
and from this point onward, data and MC calculations 
are treated identically by the analysis programs. 

All of the major ^-induced backgrounds are con- 
strained by our measurements outside the signal region. 
The inclusive CC background is verified by comparing 
data to MC calculations for events with two subevents, 
where the second subevent has < 200 tank hits and is 
consistent with a muon-decay electron. As the probabil- 
ity for fi~ capture in the oil is 8%, there are an order 
of magnitude more CC inclusive scattering events with 
two subevents than with only one subevent, so that this 
background is well checked. These data events are also 
modified by moving the hits of the second subevent ear- 
lier in time to model early, inseparable decays which can 
look more like an electron. 

To determine the NC ir° background, 7r° rates are mea- 
sured in bins of momentum by counting events in the 77 
mass peak. The MC simulation is used to correct the 
production rate for inefficiency, background and resolu- 
tion (corrections are ~ 10%). To match the data angular 
distribution, the 7r° candidates are fit to MC templates 
(in mass and angle) for resonant and coherent production 
(generated using the model of Rein and Sehgal as 
well as a template for non-7r° background events. The 
fitted parameters are used to reweight w° from the MC 
calculations and to constrain the A — > N-f rate, which 
has a branching ratio at the peak of the A resonance of 
0.56%. NC coherent 7 background and NC radiative 
7 background [3l[ are both estimated to be negligible. 
The background from interactions in the dirt surround- 
ing the detector is measured from a sample of inward- 
pointing events inside the tank at high radius. 

A sample of ^10 5 candidate z/ M CCQE events is ob- 
tained by requiring a /x-decay electron with a recon- 
structed vertex consistent with the estimated endpoint of 
the parent muon's track (60% efficiency). The observed 



rate of these CCQE events is used to correct the MC 
predictions for v e signal events, CC backgrounds, and 
v e from /1 backgrounds (which share their 7r parentage 
with the Vfj, CCQE events). These constraints increase 
the event normalization by 32% and greatly reduce the 
rate uncertainties on these three components of the final 
analysis sample. 

Systematic errors are associated with neutrino fluxes, 
the detector model, and neutrino cross sections. The 
neutrino flux systematic errors are determined from the 
uncertainties of particle production measurements, the 
detector model systematic errors are mostly determined 
from fits to MiniBooNE data, and the neutrino cross sec- 
tion systematic errors are determined from MiniBooNE 
data as well as from external sources, both experimen- 
tal and theoretical. These groups of errors are taken to 
be independent, and, for each, an individual error ma- 
trix is formed that includes the full correlation among 
the systematic parameters. This is mapped to a matrix 
describing the correlated errors in predicted background 
plus possible signal in eight v e E® E bins. The final co- 
variance matrix for all sources of uncertainty (statistical 
and systematic) is the sum of the individual error ma- 
trices. The signal extraction is performed by computing 
the x 2 comparing data to predicted background plus a 
(sin 2 (28), Am 2 )-detcrmincd contribution from — > v e 
two-neutrino oscillations in the eight E® E bins and min- 
imizing with respect to these two oscillation parameters 
across their physical range. 

With the analysis cuts set, a signal-blind test of data- 
MC agreement in the signal region was performed. The 
full two-neutrino oscillation fit was done in the range 
300 < E$ E < 3000 MeV and, with no information on 
the fit parameters revealed, the sum of predicted back- 
ground and simulated best-fit signal was compared to 
data in several variables, returning only the % 2 . While 
agreement was good in most of the comparisons, the E v i a 
spectrum had a x 2 probability of only 1%. This triggered 
further investigation of the backgrounds, focusing on the 
lowest energies where ^-induced backgrounds, some of 
which arc difficult to model, arc large. As part of this 
study, one more piece of information from the signal re- 
gion was released: unsigned bin-by-bin fractional discrep- 
ancies in the E V i S spectrum. While ambiguous, these re- 
inforced suspicions about the low-energy region. Though 
we found no specific problems with the background es- 
timates, it was found that raising the minimum E$ E of 
the fit region to 475 MeV greatly reduced a number of 
backgrounds with little impact on the fit's sensitivity to 
oscillations. We thus performed our oscillation fits in the 
energy rang e 475 < E$ E < 3000 MeV and opened the 
full data set. 

The top plot of Fig. [2] shows candidate v e events as a 
function of E® E . The vertical dashed line indicates the 
minimum Ef? E used in the two-neutrino oscillation anal- 
ysis. There is no significant excess of events (22 ± 19 ±35 
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events) for 475 < E® E < 1250 MeV; however, an ex- 
cess of events (96 ± 17 ± 20 events) is observed below 475 
MeV. This low-energy excess cannot be explained by a 
two-neutrino oscillation model, and its source is under 
investigation. The dashed histogram in Fig. [5] shows the 
predicted spectrum when the best-fit two-neutrino oscil- 
lation signal is added to the predicted background. The 
bottom panel of the figure shows background-subtracted 
data with the best-fit two-neutrino oscillation and two 
oscillation points from the favored LSND region. The 
oscillation fit in the 475 < E® E < 3000 MeV energy 
range yields a \ 2 probability of 93% for the null hypoth- 
esis, and a probability of 99% for the (sin 2 26 = 10 -3 , 
Am 2 = 4 eV 2 ) best- fit point. 
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FIG. 2: The top plot shows the number of candidate v e events 
as a function of E® E . The points represent the data with sta- 
tistical error, while the histogram is the expected background 
with systematic errors from all sources. The vertical dashed 
line indicates the threshold used in the two-neutrino oscilla- 
tion analysis. Also shown are the best-fit oscillation spec- 
trum (dashed histogram) and the background contributions 
from Vfj, and events. The bottom plot shows the number of 
events with the predicted background subtracted as a func- 
tion of E® E , where the points represent the data with total 
errors and the two histograms correspond to LSND solutions 
at high and low Am 2 . 

A single-sided raster scan to a two neutrino 
appearance-only oscillation model is used in the energy 
range 475 < E$ E < 3000 MeV to find the 90% CL limit 
corresponding to A x 2 = x\ m t ~ Xbestfit = 1M - As 
shown by the top plot in Fig. H the LSND 90% CL al- 
lowed region is excluded at the 90% CL. A joint analysis 
as a function of Am 2 , using a combined \ 2 of the best 
fit values and errors for LSND and MiniBooNE. excludes 



at 98% CL two-neutrino appearance oscillations as an 
explanation of the LSND anomaly. The bottom plot of 
Fig. shows limits from the KARMEN 0] and Bugey 
[32J experiments. 

A second analysis developed simultaneously and with 
the same blindness criteria used a different set of recon- 
struction programs, PID algorithms, and fitting and nor- 
malization processes. The reconstruction used a simpler 
model of light emission and propagation. The PID used 
172 quantities such as charge and time likelihoods in an- 
gular bins, M 77 , and likelihood ratios (electron/pion and 
electron Vmuon) as inputs to boosted decision tree algo- 
rithms [33J that are trained on sets of simulated signal 
events and background events with a cascade-training 
technique (3^ |. In order to achieve the maximum sen- 
sitivity to oscillations, the j^-CCQE data sample with 
two subevents were fit simultaneously with the t/ e -CCQE 
candidate sample with one subevent. By forming a \ 2 
using both data sets and using the corresponding covari- 
ance matrix to relate the contents of the bins of the two 
distributions, the errors in the oscillation parameters that 
best describe the z^-CCQE candidate data set were well 
constrained by the observed z^-CCQE data. This pro- 
cedure is partially equivalent to doing a v e to ratio 
analysis where many of the systematic uncertainties can- 
cel. 

The two analyses are very complementary, with the 
second having a better signal-to-background ratio, but 
the first having less sensitivity to systematic errors from 
detector properties. These different strengths resulted in 
very similar oscillation sensitivities and, when unblindcd, 
yielded the expected overlap of events and very similar 
oscillation fit results. The second analysis also sees more 
events than expected at low energy, but with less signif- 
icance. Based on the predicted sensitivities before un- 
blinding, we decided to present the first analysis as our 
oscillation result, with the second as a powerful cross- 
check. 

In summary, while there is a presently unexplained 
discrepancy with data lying above background at low 
energy, there is excellent agreement between data and 
prediction in the oscillation analysis region. If the oscil- 
lations of neutrinos and antineutrinos are the same, this 
result excludes two neutrino appearance-only oscillations 
as an explanation of the LSND anomaly at 98% CL. 
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FIG. 3: The top plot shows the MiniBooNE 90% CL limit 
(thick solid curve) and sensitivity (dashed curve) for events 
with 475 < E® E < 3000 MeV within a two neutrino oscilla- 
tion model. Also shown is the limit from the boosted decision 
tree analysis (thin solid curve) for events with 300 < E® E < 
3000 MeV. The bottom plot shows the limits from the KAR- 
MEN 0] and Bugey [H experiments. The MiniBooNE and 
Bugey curves are 1-sided upper limits on sin 2 26 correspond- 
ing to Ax 2 = 1.64, while the KARMEN curve is a "unified 
approach" 2D contour. The shaded areas show the 90% and 
99% CL allowed regions from the LSND experiment. 



ware in the analysis of the data. 
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